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Abstract 
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1 Introduction 



Extended stable or long-lived objects like kinks, vortices or domain walls are 
interesting for several reasons. First, they are copiously produced in con- 
densed matter systems, especially in time dependent conditions. Second, 
they have non-perturbative origin, often related to existence of topologically 
nontrivial sectors in pertinent models. Resulting complexity of their static 
characteristics, of evolution and interactions poses a challenge to theorists as 
well as to experimenters. Third, they likely were produced in phase transi- 
tions in the early Universe and contributed to later evolution of it. Finally, 
physics of the extended objects inspires avant-garde ideas, for example, about 
the four-dimensional Universe being a domain wall. All these aspects are 
discussed in numerous papers. Introductory lectures and extensive lists of 
references can be found in, e.g., 0, [|, f|, Q. 

In the present paper we concentrate on topological kinks. In a typical 
case one describes them in terms of a field theoretical model which has a non 
connected vacuum manifold. Basic kink is static. The pertinent solutions of 
field equations smoothly interpolate between a pair of classical vacua belong- 
ing to different connected components of the vacuum manifold. Typically, 
the fields approach the classical vacua exponentially. From the static kink 
one can get a moving one by Lorentzian or Galilean boosts. 

The topological kinks are particularly important because they are related 
to domain walls. Consider a stable planar domain wall, homogeneous in the 
two directions along the domain wall. The corresponding fields depend only 
on a coordinate £ parametrizing the direction perpendicular to the domain 
wall. They are precisely the same as for the kink. Thus, the kink can 
be identified with the transverse profile of the planar domain wall. The 
planar domain wall is the main prerequisite of analytical description of curved 
domain walls with a help of expansion in width, see, e.g., 0. 

We consider a field theoretical model which involves a single, real scalar 
field with double degerate classical ground state ±0 O - This is a typical set 
up for theoretical investigations of domain walls. However, in our case the 
field potential U(0) is not smooth. It is merely continuous, and it steeply 
increases to infinity at the ground state values ±</> of the field. At first 
glance such a model might look quite artificial. As a matter of fact, it is just 
the opposite — such field potential has very simple mechanical realization 
described in next Section. One can easily perform experiments in which the 
one-dimensional kinks are observed. This is first attractive feature of the 
presented model. 

Another interesting fact has to do with the theoretical side of the model. 
It turns out that the kink in our model has the peculiar feature that at 
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certain finite distance £o from the center of it the scalar field reaches exactly 
the classical ground states ±0 O and remains constant at larger distances. 
The common exponential approach to the ground state is absent here. 

Solutions of non-linear evolution equations which have compact support 
are called compactons. They were obtained for various modified Korteweg- 
de Vries equations |7|, §, |], [T(| pT| . The compacton presented in our paper 
is distinguished by two features. First, it has a nonvanishing topological 
charge related to the degenerate classical ground state, and therefore its 
embedding in a higher dimensional space gives a stable domain wall. In 
the K-dV case, higher dimensional generalizations are not straightrorward. 
Second, our compacton can be observed in the simple mechanical system. 

The plan of our paper is as follows. In Section 2 we discuss the topological 
compacton. The mechanical system with compactons is presented in Section 
3. Section 4 contains several remarks. 

2 The topological compacton 

We consider one dimensional topological kinks. The space coordinate is de- 
noted by £, and the time coordinate by r. Moreover, r and £ are dimen- 
sionless. The reason for this notation is that in theoretical analysis of the 
mechanical model presented in next Section r and £ appear as rescaled ordi- 
nary time t and position x, respectively. 
Lagrangian of our model has the form 

^ = ^(^) 2 -^(^) 2 -W), (1) 

where V((p) is the field potential. In order to keep the following discussion 
simple, we assume that the ground state in the model is double degenerate, 
and that it is obtained for = ±0 O . The field potential is symmetric, 
V((p) = V(—(j)), and it has a local maximum at = 0. Considerations 
presented below can easily be generalized to fields potentials of other forms. 

As always, the kinks are given by solutions of the field equation obtained 
from Lagrangian (1). They interpolate between the two ground state values 
of the field. For the static kink the field equation has the form 

8$<f> - V\<P) = (2) 

where V = dV/dcf). Multiplying Eq.(2) by and integrating over £ we 
obtain the equation 

i(<%0) 2 - V(<j>) = constant. (3) 
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We expect that d^(f> — > when — * 0o- Therefore, constant = —V((j)o), and 
Eq.(3) can be written in the form 

(^0) 2 = 2(\/(0)-\/(0 o )). (4) 

For the kink, 0(0 increases with £ monotonically from — 0o to +0o- Hence, 
c?£0 > 0, and Eq.(4) can be written as 

<9 5 = p (V(<j>)-V(<j>o)). (5) 

This equation is very convenient starting point for a detailed analysis of the 
kink. 

Let us first have a look at the behaviour of when it approaches 0o- 
Then, 

where <50 > 0. Expanding the l.h.s. of Eq.(5) in S(j) and keeping only the 
leading term we obtain the following equation 



2|\/'(0 O )|V50. (6) 

Here V\4>o) is °f course understood as the limit from the side of < O (the 
same applies to all derivatives V^(0 O ) below). General solution of Eq.(6) 
has the form 

50(o = ^in0o)i(£ O -o 2 , (7) 

where £o is an arbitrary constant. Thus, we have obtained the parabolic 
approach to the ground state value of the field O . This value is reached at 
£ = £ exactly. 

The unusual polynomial approach of the field to its ground state value is 
of course due to the fact that V'((f>o) ^ 0. Then, in the expansion 

V(0)-V(0 o ) = n0o)(0-0o) ( 8 ) 
+ V(0 O )(0 - 0o) 2 + ^'"(0 O )(0 - O ) 3 + • • • 

the first term dominates the limit — > 0o~ . The well-known exponential 
approach is obtained when V'((p ) = and V"((f> ) > 0. In this case 



50(O = exp(-^"(0 o )O. (9) 

If the first nonvanishing term is of the order k > 2 and (-l) fc \/( fc )(0 o ) > 
then the approach to 0o is power-like with a negative power of £, namely 



50(0 



k / 2|\/W(0q)| 
2 _1) V k\ 



2 

2-k 



2 

£ fc -2 



(10) 
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The power-like asymptotics (10) is sufficient for convergence at large £ of the 
integral in the expression for the total energy E k of the kink, 



/+oo / 1 



(ii) 



The term V((f>o) is included in order to subtract the energy of the ground 
state. 

The kink located at £ = is the solution 0(£) of Eq.(5) which has the 
following dependence on £: 0(£) is equal to — 0o for £ < — £ , then it starts 
to increase in the parabolic manner, 

0(O = -0o + ^IH0o)|(£ + £o) 2 , 

it reaches at £ = 0, and continues to increase to +0o as £ approaches No- 
dose to £o 

0(O = 0o-^in0o)|(£-£o) 2 . 

Finally, 0(£) = 0o for £ > £o- Anti-kink is obtained by changing the sign of 
0(£). It is clear that such kinks and anti-kinks are compactons. 
As the example, let us take the field potential V(4>) of the form 



v(4>) = 



COS ( 

oo 



1 for 
for 



< 
> 



(12) 



where 0o is a constant, n > O > 0, see Fig.l. It has two degenerate minima 
at 4> — =t0o- 




Fig.l. Plot of the potential V(4>) for O — 1-1- 
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Because of the infinite potential barrier, modulus of the field cannot 
exceed O . Therefore, Eq.(5) is now supplemented by the condition 



~ 00 < < 00- (13) 

Potential of this kind has the simple mechanical realisation described in next 
Section. The angle 0o is an external parameter which we can easily regulate. 

The length of the compacton is equal to 2£ - It is related to 0o by the 
following formula obtained by integration of Eq.(5) 



r4>o 

/ d<P- 
Jo 



\J COS — COS 0o 

The change of the integration variable to A = 0/0 o gives formula 



Jo 



'/A , / \ = ^o, (14) 

cos(A0 o ) — cos 0o 



o we may expand the cosine functions. Keeping the terms up to ~ 0q we 



from which we can find the dependence of £o on O . In particular, for small 
0o we may expand the cosin 
obtain the following formula 

6 = ^(1 + ^ + ...). (15) 




Thus, the length of the kink remains finite in the limit 0o — > 0. Nevertheless, 
in this limit the kink of course disappears because it becomes completely flat, 

0(O~o. 

Because 0(— £) = —0(0: it i s sufficient to obtain the kink solution in the 
half-line £ > 0. The initial value of at £ = is known: 0(0) = 0. Then, 
Eq.(5) is equivalent to the following integral equation 



L 



C (16) 



ycos(0 o A) — cos 0o 



The integral on the l.h.s. can be related to an elliptic function. The integral 
equation (16) can be used for calculation of the approximate forms of 0(0- 
For small £ we just expand in A the cosine function cos(0oA) and the inverse 
of the square root. We find that for £ — > 

0(O = 2sin^-isin^ + .... (17) 
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In order to obtain 0(£) for £ — > £o _ we subtract from both sides of Eq.(16) 
a/2£o/0Oi an d on the l.h.s. we use formula (14). We obtain the integral over 
A in the interval [0(0/00) !]■ Such integral can be tackled by expanding the 
integrand in e = 1 — A. Finally, we find that for £ — > £ _ 

^ a sin0o . v 2 sin(20 o ) 4 

0(O = 0o (f-fo) ^ — (f-fo) +•••• (18) 

Also obtaining a numerical solution is rather easy. The solution 0(£) 
starts from = at £ = with the slope 



d, 



. 00 



'tv\z=o = 2sin y' 



and it has the smooth shape presented in Fig. 2. It has been obtained with 
the help of Maple©. 




Fig. 2. The numerical solution 0(£) of Eq.(5). V(<f>) is given by formula (12) 
with 0o = 1.1. 



The total energy of the field is given by the formula 
Using Eq.(4) we can write it also as 

-a 
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e = r° ^^(^0) 2 +cos0-i 

J — f n L Z 



E = -2f (l - cos0 o ) + / & d£(%0) 2 . 



(19) 



(20) 



The first term on the r.h.s. of formula (2) gives just the energy of the back- 
ground on the segment occupied by the kink. The second term is equal to 
the proper energy of the kink. It can be regarded as its rest mass M. We 
have computed E and M numerically. Sample results are given in the Table 
below. 



00 


E 


M 




1.5 


-0.263 


3.128 


1.824 


1.0 


-0.062 


1.478 


1.674 


0.2 


-0.019 


0.044 


1.574 


0.05 


-0.002 


0.002 


1.570 



Table 1. Results of numerical computations of the total energy of the field 
and of the rest mass of the compacton. 

We see that the total energy E of the field is negative in spite of the 
presence of the kink with the positive rest mass. We have checked that this 
is the case also for other choices of 4>o- 

3 The mechanical system with the topological 
compacton 

Let us take a thick rectilinear wire with 2N + 1 pendulums connected to it 
at the points x iy % = —N, ...,N. The points Xj are separated by a constant 
distance a. Each pendulum has a very light arm of length R and a mass m 
at the free end. It can swing only in the plane perpendicular to the wire. All 
pendulums are fastened to the wire stiffly, hence their swings twist the wire 
accordingly. The wire is elastic with respect to such twists. Each pendulum 
has one degree of freedom which may be represented by the angle 
between the vertical direction and the pendulum. Thus, $(xj,t) = corre- 
sponds to the upward position of the i-th pendulum. When all pendulums 
point upwards the wire is not twisted. 

Equation of motion for each pendulum, except for the first and the last 
ones, has the following form 

2 <P$( Xi ,t) . $(x i -a,t) + $(x i + a,t)-2$(x i ,t) 
mR — = mgRsm&(xi,t) + k , 

(Xv (X 

(21) 

where k is a constant which characterizes the elasticity of the wire with 
respect to the twisting. The first term on the r.h.s. of Eq.(21) is due to 
the gravitational force acting on the mass m, the second term is the elastic 
torque due to the torsion of the wire. 
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As for the two outermost pendulums, we assume that they are kept in 
the upward position by some external force, that is that 

$(x_at, t) = 0, $(x N ,t) = 0. (22) 

The mechanical system described above is esssentially identical with one 
used in a realization of sinus-Gordon solitons, except that in the case of 
these solitons all pendulums initially point downwards, so that they are in 
the stable equilibrium position. In the present case they are put into the 
seemingly unstable equilibrium position &(xi,t = 0) = 0. However, as we 
show below, this position can actually be stable due to the presence of the 
external force which implements conditions (22). Another difference with the 
sinus-Gordon case is that we mechanically restrict the range of <3>(xj,t) 

\$(xi,t)\<$ <ir, (23) 

by putting on both sides of the wire and parallely to it two rigid rods. The 
pendulums rebound from the rods if $(xj, t) = ±$o- 

Let us introduce a function $(x, t), where a; is a real continuous variable, 
twice differentiable with respect to x and such that its values at the points 
equal to the angles $(xi,t). Hence, $(x, t) interpolates between 
$(xj, t). Of course, for a given set of values of the angles there is an infinite 
number of such functions. For any of them the following identity holds 



ra rO 

®(xi — a,t) + $(xi + a,t)—2$(xi,t)= / ds 1 / ds 

JO J -a 



<9 2 $(S! + S 2 +X,t) 
2 



dx 2 



We shall restrict our considerations to such motions of the pendulums that 
there exists the interpolating function $(x, t) of continuous variables x, t such 
that 

2 d 2 Hx,t) 



/ d Sl / 

JO J-a 



o J <9 2 $( Sl + s 2 + x,t) 
ds 2 



dx 2 



a 



dx 2 



(24) 



for all times t and at all points Xi. For example, this is the case when the 
second derivative of $ with respect to x is almost constant when x runs 
through the interval [x« — a, Xi + a], for all times t. With the approximation 
(24) the identity written above can be replaced by the following approximate 
one 

$(xi -a,t) + $(xi + a,t)- 2$(xi, t) w a 2 9 ^ 

Using this formula in Eq.(21) we obtain 

~d 2 $(xi,t) „ _ . d 2 $(x,t) 

mR 2 v 1 ~ mgR sin $(^ , t) + kcl 



dt 2 * v " ' dx 2 



(25) 
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Let us now suppose that our function t) obeys the following partial 
differential equation, 

,<9 2 $(a;,t) „ _ . <9 2 $(:r,t) , x 

mR 2 = mgRsm$(x,t) + kcl — tt^, (26) 

at 2 ox 2 

where x G [—Na,Na], and 

(f)(-Na,t) = 0, (f)(Na,t) = 0. (27) 

Then, it is clear that <3>(xj,t) obey Eq.(25) and the boundary conditions (22). 
Hence, if condition (24) is satisfied we obtain the approximate solution of the 
initial Newton equations (21). In the final step, we pass to the dimensionless 
variables 

V K \ KCL 

Then, Eq.(26) acquires the form 

d 2 0(£,r) d 2 0(£,r) 



dr 2 oe 



sin0(£,r) = O. (28) 



This equation coincides with the Euler-Lagrange equation obtained from La- 
grangian (1) with potential (12) for |0| < O - The restriction (23) now has 
the form 

I0(£,r)|<0 o . 

The range of the spatial coordinate £ is from — £tv to £at, where £at = 
N^JingRa/n. The boundary conditions for have the form 

0(-&v,r) = O = 0(£v,t). (29) 

Equation (28) with boundary conditions (29) has the trivial solution 

0(£,r)=O. (30) 

Let us check stability of this solution against small perturbations. To this 
end we write <f>(£, r) = e(£, r) and we linearize Eq.(28) with respect to e. The 
resulting equation has the form 

It leads to the following dispersion relation 



7T 



u(k) = ±\Jk 2 — 1, k = -^—n, 
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where n is an integer different from zero. The discrete values of k are due 
to the boundary conditions (29). We see that the trivial solution is stable 
if £jv < 7r/2. Loosely speaking, this condition is satisfied when the wire is 
short enough. 

The kink can appear only if the trivial solution is not stable. Then the 
boundary conditions (29) can be satisfied by locating one half of the kink at 
the one end of the wire and one half of the anti-kink at the other end. Each 
of them has the length £ , where £o has been introduced in previous Section. 
We have to put at least one full kink in the middle. Therefore, the total 
dimensionless length of the wire has to be equal at least to 4£ . Thus, full 
compactons can be observed if 

6v > 2£o. 

Finally, let us note that we made a mechanical model of that kind, and 
the compactons were clearly seen. 

4 Remarks 

1. The single kink solution discussed in Section 2 can be trivially general- 
ized. Lagrangian (1) and condition (13) are Poincare invariant. Therefore, 
Lorentzian boosts give kinks moving with arbitrary velocities not exceeding 
1 in the dimensionless variables. Because the kinks do not feel each other 
when their centers are separated by a distance larger that 2£ , it is clear that 
there exist static multi-kink solutions. Such 'trains' of separated kinks and 
anti-kinks are also the compacton type solutions. Moreover, also kinks and 
anti-kinks moving with different velocities can be trivially combined to give 
compacton solutions in a finite time interval until they touch each other. 

2. Our work can be extended in several directions. One could study the 
related non-planar domain walls. We have mentioned in the Introduction 
that they can be constructed with the help of expansion in width. 

Another direction is suggested by the fact that the separated compactons 
do not interact. Therefore, they seem to be well-suited for testing theoret- 
ical descriptions of productions of topological defects, 0, |], |}]. The lack of 
interactions in the final state simplifies counting of the defects. Our current 
work is devoted to this topic. 

One could also study interaction of the kink with the anti-kink, for ex- 
ample, when they scatter on each other with various relative velocities. 



11 



5 Acknowledgement 

We would like to thank Dr. P. W§grzyn for remarks pertinent to this work. 

References 

[1] P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter 
Physics. Cambridge University Press, 1995. 

[2] T. W. B. Kibble and M. B. Hindmarsh, Rep. Prog. Phys. 58, 447 (1995). 

[3] Y. M. Bunkov and H. Godfrin (Eds.), Topological Defects and the Non- 
Equilibrium Dynamics of Symmetry Breaking Phase Transitions. Kluwer 
Academic Publishers, 2000. 

[4] W. H. Zurek, Phys. Rep. 276, 177 (1996). 

[5] M. Pavsic, The Landscape of Theoretical Physics: a Global View. From 
Point Particles to the Brane World and Beyond, in Search of a Unifying 
Principle. Kluwer Academic Publishers, 2001. 

[6] H. Arodz, Acta Phys. Pol. B29, 3725 (1998). 

[7] P. Rosenau and J. M. Hyman, Phys. Rev. Lett. 70, 564 (1993). 

[8] F. Cooper, H. Shepard and P. Sodano, Phys. Rev. E48, 4027 (1993). 

[9] A. Khare and F. Cooper, Phys. Rev. E48, 4843 (1993). 

[10] B. Dey and A. Khare, Phys. Rev. E58, R2741 (1998). 

[11] F. Cooper, J. M. Hyman and A. Khare, Phys. Rev. E64(2) , 026608 
(2001). 



12 



